arXiv:astro-ph/0503067vl 3 Mar 2005 


Astronomy & Astrophysics manuscript no. 2227 February 2, 2008 

(DOI: will be inserted by hand later) 


A Field Programmable Gate Array Spectrometer for Radio 

Astronomy 

First Light at the Effeisberg 100-m Teiescope 

S. Stanko', B. Klein^, and J. Kerp' 

* Radioastronomisches Institut der Universitat Bonn, Auf dem Htigel 71, D-53121 Bonn 
e-mail: stanko@astro.uni-borm.de e-mail: jkerp@astro.uni-bonn.de 
^ Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn 
e-mail: bklein@inpifr-boim.mpg.de 

Received 22/10/2004; accepted 22/02/2005 

Abstract. We describe the technological concept and the first-light results of a 1024-channel spectrometer based on field 
programmable gate array (FPGA) hardware. This spectrometer is the prototype for the seven beam L-band receiver to be 
installed at the Effeisberg 100-m telescope in autumn 2005. Using “of-the-shelf” hardware and software products, we designed 
and constructed an extremely flexible Fast-Fourier-Transform (FFT) spectrometer with unprecedented sensitivity and dynamic 
range, which can be considered prototypical for spectrometer development in future radio astronomy. 
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1. Introduction 

A wealth of information on the physical conditions of objects 
in the universe can be gathered using spectrometers. In the ra¬ 
dio regime we can roughly identify three basic types of spec¬ 
trometers: autocorrelators, acusto-optical spectrometers and fil¬ 
ter banks. Today these spectrometers offer a useable bandwidth 
from a few kHz up to 2 GHz with a few thousand spectral chan¬ 
nels, which are capable of resolving narrow spectral lines of 
masers and the thermal line emission of gaseous clouds. 

Over the past decade radio astronomy has been in a tran¬ 
sitional phase from single to multifeed arrays in continuum 
(bolometers) as well as in heterodyne receivers (spectroscopy). 
Multifeed receivers enable very sufficient surveys to large sky 
areas and significantly improve the signal-to-noise ratio using 
the on-the-fly observing mode. The huge impact of a multifeed 
receiver on H i and pulsar astronomy (Manchester et al. I200H 
can be well illustrated by the Parkes 21-cm multifeed system 
(Stavelev-Smith n~996t . Unfortunately, the backends do not de¬ 
velop at the same pace as the receiver technology. Neither the 
number of spectral channels nor the number of spectrometers 
nor the bandwidths have increased significantly. 

FPGA-chips might offer a solution for this “backend bottle 
neck”. Using megagate chips, FFT- or polyphase filter bank al¬ 
gorithms can be implemented today with spectral channel num¬ 
bers ranging between 16,000 to 32,000. The huge advantages 
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of FPGAs are low price, wide availablity (both based on the 
large commercial interest) and the “plug-and-play” installation 
as commercial personal computer (PC) - PCIbus cards, which 
offer a high reliability. Moreover, the sub-mm observatories at 
high altitudes have severe difficulties to provide the receivers 
and supplementary power aggregates, especially for digital au¬ 
tocorrelators that need several kW. FPGA-spectrometers con¬ 
sume only a few W, so heat transfer via coolers is unneces¬ 
sary. The low power consumption and compact size of FPGA- 
spectrometers are attributes that recommends this type of back¬ 
end for use on spacecraft and satellites. 

In this paper we present our design approach toward devel¬ 
oping an FPGA-spectrometer for radio astronomy. The tech¬ 
nical concept is presented in Sect. 2 of this research note, fol¬ 
lowed by the basic implementation. First light observations of 
the Hi and OH maser line emission at the Effeisberg 100-m 
Telescope are presented in Sect. 3. 

2. Techniques 

2.1. Concept 

One of the most important fundamentals for correctly sampling 
analog signals is the concept of aliasing. The famous Nyquist 
criterion states that data recorded at a rate of samples per 
second can effectively represent a signal of bandwidth up to 
Fs/2 Hz. Sampling signals with greater bandwidth produces 
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FFT - Spectrometer 



Fig. 1. Sketch of the ICS 554C device and the arrangement of Intellectual Property (IP)-cores to constitute an FFT-spectrometer 
on the field programmable gate array (FPGA). The two signals from the 21-cm receiver enter the device on the left-hand side 
via the analog to digital converters (ADCs). The optional digital down converters (DDCs) reduce the bandwidth to values from 
10 MHz to a few kHz. Within the FPGA, we implement a distributed halfband filter (DHBF), a window-function (look-up table) 
followed by the FFT. Both the DDC and the DHBF produce complex data values (I: in-phase, Q: quadrature). The pipeline FFT 
output is transferred into the power-spectrum calculator which sums the spectra in an averager. The spectra are stored on the PC 
via the first-in first-out pipelines (FIFOs). 



Fig. 2. Principle of downconverting in Full Band Mode and 
downconverting with additional filtering, mixing and decima¬ 
tion in Narrow Band Mode. The output frequency of the re¬ 
ceiver IF band is 150 to 200 MHz (a). Because of the sampling 
frequency of 100 MHz the 4th Nyquist Band is sampled (b). In 
DHBF mode the whole band with center frequency of 25 MHz 
ist converted to OMHz (c). In Narrow Band mode (using the 
DDCs) a part of the band (dashed in the example) with a user 
selectable center frequency is converted to 0 MHz. 


aliasing; signal components at frequencies greater than Fsl2 
are folded, or aliased, back into the FJl band. 

The frequency translation aspect of aliasing can be ex¬ 
ploited by using a technique called undersampling, band-pass 


sampling, harmonic sampling or Super-Nyquist sampling. To 
understand this technique, one must carefully consider the def¬ 
inition of the Nyquist constraint. Sampling a signal of band¬ 
width Fj/2 requires a minimum sample rate Fs, but this FJl 
bandwidth can theoretically be located anywhere in the fre¬ 
quency spectrum (e.g. nF^ to (n -i- 1/2) n - 0,1,2,...), 
not simply from DC to FJ2 like baseband signals. The alias¬ 
ing action, like a digital mixer, can be used to translate an in¬ 
termediate frequency (IF) down to the baseband (see Fig.|2}. 
Signals s in the bands nFs < s < {n 1/2) Fj will be 
translated down intact; the translation of signals in the bands 
(n - 1 /2) Fs < s < nFs will be “mirrored” in frequency, an 
effect corrected in the PC after the measurement. 

Using undersampling techniques (Groshong & Ruscak 
I1991t enables direkt IF-band sampling of the receiver outputs 
without any additional analog baseband mixing. Only one (an¬ 
tialiasing) bandpass filter in front of the analog-digital con¬ 
verter (ADC) is necessary to select the Nyquist band and avoids 
aliasing from other frequency bands. This concept saves costs 
and, because analog mixers are highly sensitive to input levels, 
increases the signal quality. If not carefully adjusted, the input 
signal can cause considerable harmonic distortion, which, in 
addition varies as a function of temperature. 

The system is configured to process two channels in par¬ 
allel, each channel with a 50 MHz bandwidth. It is based on a 
data acquisition card ICS554C (ICS I2003t . consisting of two 
14-bit ADCs running at 100 MHz sample rate, a Xilinx FPGA, 
two FIFO-buffers and a PCIbus interface (Fig.^l- 

The on-board Xilinx XC2V3000 FGPA is capable of pro¬ 
cessing the signal from both ADCs in parallel to produce 
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2 X 1024 channel spectra. The digital data streams produced by 
both ADCs are connected to the FPGA and piped in through 
separated Digital Down Converters (DDC) chips which are in¬ 
stalled to determine the bandwidth (increase the spectral res¬ 
olution) of the signal and to select sub-bands. Any necessary 
mixing and filtering is done digitally in these DDCs. No ana¬ 
log elements are used and reconfigurating the sub-band and 
bandwidth is easily done via the PC. The DDCs are GC4016 
quad receivers (Gravchip l2000> . each one containing four iden¬ 
tical down conversion circuits. The input data stream can be 
received in four separate sub-bands, each with up to 2.5 MHz 
bandwidth. By combining all of these, a maximum bandwidth 
of 10 MHz can be processed, corresponding to a spectral reso¬ 
lution of about 10 kHz per spectral channel. The best resolution 
of ~ 20 Hz per channel can be achieved by reducing the mini¬ 
mum bandwidth to 20kHz. 

2.2. Implementation 

Individual cores were purchased from RF-Engines to imple¬ 
ment the processing pipeline shown in Fig.[2 This straightfor¬ 
ward implementation and concept makes it possible to obtain 
debugging information at any stage of the pipeline. 

The “raw” data from the ADCs, as well as the data pro¬ 
cessed by the DDCs, are routed into the FPGA in parallel. A 
virtual switch has been implemented in order to select data 
either as Full Band (50 MHz bandwidth) from ADC (also re¬ 
ferred to as the Distributed Halfband Filter, DHBF) or Narrow 
Band (10 MHz - 20kHz bandwidth) from the DDCs. When 
using the DHBF, the bit stream is mixed with sine and cosine 
waves to shift the signal, centered symmetrically about the IF 
of half the bandwidth, to a signal centered at 0 Hz (see Fig.|2] 
(b-c)). In contrast to real-only digitization this technique yields 
a full Nyquist sampling, and the following pipeline FFT works 
with the highest efficiency on the complex data. 

An additional FIFO had to be added to deal with data by 
the DDCs. This is because the pipeline always transfers data of 
exactly 2048 samples, but the data produced by the DDCs are 
not delivered continuously. 

A window function is applied to the data first (Kaiser n'974t . 
This function is freely programmable via the PC without re¬ 
loading the FPGA design by means of “slow” PCIbus connec¬ 
tion. 

The central element of the whole processing pipeline is the 
Vectis HiSpeed 1024-point radix-2 FFT algorithm developed 
by RF-Engines. This pipelined EET-core works up to 100 MHz, 
enabling a continuous Eourier Transformation of the input sig¬ 
nals. The output consists of 27-bit interleaved real and imagi¬ 
nary signed data, which are fed into the power spectrum calcu¬ 
lator and the following averager. 

There are two possible communication channels between 
the EPGA and the PC: on the one hand a “slow” 8 bit PCI con¬ 
trol bus, on the other hand the ICS554C card with two outgoing 
hardware EIEOs, which can transfer data via a 64-bit, 64 MHz 
PCI bus. This gives a maximum PCI burst data transfer speed 
ofupto528 MBs-*. 
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Fig. 3. Allan-plot of the stability of the EPGA-spectrometer. 



time [s] 

Fig. 4. Effelsberg 100-m RMS-limit temperature as a function 
of integration time of the 1024 and 8192 spectral channel auto¬ 
correlators with the EPGA-spectrometer. The autocorrelators 
are the standard spectrometers at the Effelsberg 100-m tele¬ 
scope. The 1024 channel autocorrelator samples 1-bit, the 8192 
channel autocorrelator 2-bit/3-level. 

The control bus is used to set window, integration time 
(number of samples) and observing mode, as well as a read¬ 
out of data ready and status information. The EIEOs are used 
to transfer the Eourier transformed and pre-integrated spectra. 
The high speed connection is essential to keep the sub-samples 
integration time low, e.g. for radio frequency interference (REI) 
investigations or pulsar measurements. 

PC software was developed that is capable of setting con¬ 
trol values for the window, integration time and the configura¬ 
tion of the DDCs. The PC is responsible for triggering the first 
sample and then wait for a data ready signal. The PC software 
then reads back the integrated spectra from 2 channels using 
the hardware EIEOs in EPGA and DMA mode. The data are 
stored as ASCITfiles with an optional header. 

2.3. Stability Tests 

In order to quantify the stability of the spectrometer, we ex¬ 
amined its performance when connected to a wideband noise 
generator. The Allan plot (Schieder 119851 generated by this 
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Fig. 5. Top: Hi 21-cm line spectrum of the standard calibra¬ 
tion source S7 observed using the total bandwidth of 50 MHz 
of the FPGA-spectrometer. The hardware bandpass was not op¬ 
timized to be flat across the whole frequency range. Bottom: 
ON- OFF subtracted DDC H i spectrum of the S7 source with 
1 minute of integration time. The bandwidth is about 2 MHz. 

exercise implies us stability times « 1000 seconds (see Fig. El- 
Compared to other spectrometers this is excellent and highly 
sufficient for astronomical applications. 

The first quantitative evaluation of the lAU standard posi¬ 
tion S7 H i-line spectra (Kalberla et al. IT^ suggested that the 
new FPGA-spectrometer using samplers with 14-bit dynamics 
is a factor of about 7r/2 more sensitive than the 1024-channel 
autocorrelator, which uses 1-bit sampling. A few days later 
in a second observing run, we were able to observe S7 with 
a cold 21-cm receiver. We compared the FPGA-spectrometer 
data with the 1024- and 8192-autocorrelator (which uses 2- 
bit/3-level sampling) spectra observed in parallel. Our results 
are plotted in Fig.|4] 

The sensitivity of the system, given by the standard ra¬ 
diometer formula 

AT C 

f sys yjB tint 

with C = I for 1-bit autocorrelators, C = 1.23 for the 2-bit/3- 
level autocorrelator, and C ^ 1 for the FPGA-spectrometer 
(Hagen & Farley I1973> . clearly indicates superior sensitivity 
of the FPGA-spectrometer. The fit gives us a ratio of 1.53 be¬ 
tween the FPGA-spectrometer and the 1024 channel / 1-bit au¬ 
tocorrelator C values and a ratio of 1.21 between the FPGA- 
spectrometer and the 8192 channel autocorrelator / 2-bit, 3- 
level. 

3. Observation 

The first light observation using our FPGA-spectrometer oc- 
cured on 24 August 2004. To study the quality of the Hi 21- 
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Fig. 6. Top: OH-maser emission of the W30H star formation 
complex observed with a bandwidth of about 2 MHz. The in¬ 
tensity contrast between the main emission line at 1667 MHz 
(right) and the satellite line at 1665 MHz (left) is an order of 
magnitude. Bottom: Resolved line emission of the OH maser 
emission at 1667 MHz (one circular polarisation channel) with 
a bandwidth of 0.2 MHz using the DDCs. Please note that the 
temperature scale is not calibrated in contrast to the H i spectra! 


cm line spectra we chose the northern sky standard calibration 
source S7. For a first look, we used the whole band (50 MHz) of 
the FPGA-spectrometer. Fig.|3(top) shows the recorded spec¬ 
trum of the standard calibration source. We intended to use 
a frequency switching mode to determine the baseline of the 
spectrometer as accurately as possible. This was done by inte¬ 
grating for 30 seconds at the center-frequency (ON spectrum), 
shifting the frequency by tuning the receiver by 6 MHz, and in¬ 
tegrating for an additional 30 seconds (OFF spectrum). To nor¬ 
malize the line emission detected by the FPGA-spectrometer, 
we used the S7-spectra observed in parallel with the autocor¬ 
relator spectrometer operated and, implementing standard H i 
data reduction procedures (Peter Kalberla, priv. comm.) de¬ 
termined the system temperature and evaluated the quality of 
the S7-calibration observation. The velocity resolution of the 
50MHzFPGA data corresponds to 10.32kms \ We computed 
the difference of the ON and OFF spectra and performed a rel¬ 
ative calibration using a part of the baseline which is free of 
astronomical line emission. 

Since the velocity channel width of 10.32 km s“' is much 
too coarse for galactic Hi studies, we switched to the DDC 
mode, setting the initial bandwidth of 2 MHz. 88% of this band¬ 
width is scientifically usable because of the digital bandpass 
of the DDCs, which suppresses the band edges. Following the 
same frequency switching procedure as described above, we 
observed the S7-line profile with a 0.5kms“^ velocity resolu¬ 
tion (Fig.|3bottom). The baseline subtraction was performed as 
described above. 
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To evaluate the dynamic-range of the FPGA-spectrometer, 
we observed the OH maser lines of the W30H-star forming 
region (Norris & Booth [TWH . Using a 2 MHz bandwidth, we 
simultaneously observed 120 K line emission at 1667 MHz and 
about an order of magnitude stronger line emission of the maser 
at 1665 MHz. Fig.|6l(top) shows these data without any base¬ 
line subtraction. All DDC spectra reveal a sinusodial baseline 
modulation, which is very stable in time. Accordingly we al¬ 
ways obtain an accurate flat baseline after a simple ON-OFF 
subtraction as described above. Using the DDCs to increase the 
velocity resolution to 0.03 km s ', we can resolve an individual 
line complex (Figl^bottom). Please note that the OH-spectra, 
in contrast to the S7 H i-spectra, are not normalized. It is worth 
noting that we do not need to change any attenuator in our setup 
when switching from the weak continuum S7 21-cm line obser¬ 
vations to the OH-maser region. This is necessary when using 
the autocorrelator spectrometer with its limited dynamic range. 


4. Conclusion 

First light observations using an FPGA-spectrometer have been 
performed at the Effelsberg 100-m Telescope. These obser¬ 
vations demonstrate that FPGA-spectrometers are capable of 
superseding curi'ent standard spectrometers with their limita¬ 
tions in bandwidth, number of spectral channels and dynamic 
range. Moreover, the sensitivity is superior to autocorrelators 
due to the larger number of ADC bits. Considering the high 
redundancy and low costs offended by commercial products, 
FPGA-spectrometers might well be the future in radio astron¬ 
omy spectroscopy. 
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